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CN ■ Abstract 

O ■ In this paper we calculate the Casimir energy for a dielectric-diamagnetic cylinder 

with the speed of light differing on the inside and outside. Although the result 



(N ■ is in general divergent, special cases are meaningful. The well-known results for a 

^ ■ uniform speed of light are reproduced. The self-stress on a purely dielectric cylinder 

P\J , is shown to vanish through second order in the deviation of the permittivity from 

I its vacuum value, in agreement with the result calculated from the sum of van der 

^ ■ Waals forces. These results are unambiguously separated from divergent terms. 
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1 Introduction 



Interest in quantum vacuum phenomena, subsumed under the general rubric 
of the Casimir effect, is increasing at a rapid pace. Status of work in the field 
is summarized in recent review articles and monographs [1,2,3]. The theoreti- 
cal developments have been largely driven by experimental and technological 
developments, where it is becoming evident that Casimir forces may present 
fundamental limits and opportunities in nanomechanical devices [4] and nano- 
electronics [5]. Thus it is imperative to understand fundamental aspects of the 
theory, such as the sign of the effect, which, at present, cannot be predicted 
without a detailed calculation. This paper represents an incremental increase 
in our list of solved examples of Casimir energies with nontrivial boundaries. 

The Casimir energy for an uniform dielectric sphere was first calculated in 
1979 by Milton [6] and later generalized to the case when both the electric 
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permittivity and the magnetic permeability are present [7]. It was later ob- 
served [8] that, in the special dilute dielectric case where fi = 1 and \s—l\ ^ 1, 
the series expansion in £ — 1 has a leading term that perfectly matches the 
"renormalized" energy obtained by summing the van der Waals interactions 
[9]. That result 

23 

is obtained either by isolating and extracting surface and volume divergences, 
or directly by analytically continuing in the number of space dimensions. 

The Casimir analysis for the case of the circular cylinder has been attempted 
on several occasions; however, the difficulty of the geometry and the fact that 
the TE and TM modes do not decouple makes the problem considerably more 
complex. Only in the case when the speeds of light inside and outside the 
cylinder coincide is the result completely unambiguous [10,11,12,13]. This in- 
cludes the classic case of a perfectly conducting cylindrical shell [14] where 
the energy per unit length is found to be 

£ = -^. (2) 

a 

where a is the radius of the cylinder. The minus sign indicates that the Casimir 
self-stress is attractive, unlike the Boyer repulsion for a sphere [15]. 

When the speed of light is different inside and outside of the body, the Casimir 
energy will be divergent [6], which goes beyond those divergences associated 
with curvature [16,17,18,19,20]. Thus it seems impossible to ascribe any signif- 
icance to results of such calculations. Any success in extracting a meaningful 
result in such cases, as in the dilute dielectric sphere example mentioned above, 
seems noteworthy. 

We present here the calculation of the Casimir pressure on the walls of an in- 
finite circular dielectric-diamagnetic cylinder with electric permittivity e and 
magnetic permeability fi inside the cylinder which is surrounded by vacuum 
with permittivity 1 and permeability 1 so e/x 7^ 1. It is shown that the cor- 
responding Casimir energy per unit length is divergent, as expected, but, for 
fjL — l, the finite coefficient of — l)^/a^ in the expansion for the dilute ap- 
proximation yields the surprising zero result found by summing the van der 
Waals energies between the molecules that make up the material, in a man- 
ner similar to that which resulted in (1) [10,21]. The latter calculation was 
independently carried out by Milonni [22], and verified by a perturbative cal- 
culation by Barton [23]. Although there should be divergences in the energy 
proportional to {e — l)^a and {e — l)^/a, the coefficient of {e — is unique 

and finite [24]. 

The paper is laid out as follows. In Sec. 2 we calculate the dyadic Green's func- 
tions that will allow us to compute the one-loop vacuum expectation values of 
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the quadratic field prodiicts. This enables us to calculate the vacuum expec- 
tation value of the stress tensor, the discontinuity of which across the surface 
gives the stress on the cylinder, computed in Sec. 3. The bulk Casimir stress, 
which would be present if either medium filled all space, is computed in Sec. 4 
and must be subtracted from the stress found in Sec. 3. Finally, the case of a 
dilute dielectric cylinder is considered in Sec. 5. and by detailed analytic and 
numerical calculations in Sec. 6 it is shown that the Casimir stress vanishes 
both in order e — 1 and — 1)^. The significance of divergences encountered 
in the calculation is discussed. The implications of these results are briefly 
considered in the Conclusions. 



2 Green's Function Derivation of the Casimir Energy 



In order to write down the Green's dyadic equations, we introduce a polariza- 
tion source P whose hnear relation with the electric field defines the Green's 
dyadic as 

E(x) = j ( dx')T{x, x') ■ P(x'). (3) 

Since the response is translationally invariant in time, we introduce the Fourier 
transform at a given frequency uj, 

T{x,x') = r ^exp[-ic.(t-0]r(r,r',c.). (4) 

We can then write the dyadic Maxwell's equations in a medium characterized 
by a dielectric constant e and a permeability both of which may be functions 
of frequency (see Ref. [6,25,26]): 



V X r'-ia;/x* = -V X 1, (5a) 
e 

-V X * -ia;£r' = 0, (5b) 
where 

r'(r,r',a;) = r(r,r',a;) + -^, (6) 
and where the unit dyadic 1 includes a three-dimensional 5 function, 

l = l(5(r-r'). (7) 

The two dyadics are solenoidal, 

V-* = 0, (8a) 

v.r' = o. (8b) 
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The corresponding second order equations are 



(V2 + a;2e/x)r' = --V X (V X 1), (9a) 

e 

(V^ + u; V)* = V X 1. (9b) 

Quantum mechanically, these Green's dyadics give the one-loop vacuum ex- 
pectation values of the product of fields at a given frequency uj, 



(E(r)E(r')) = yr(r,r'), (10a) 
(H(r)H(r')) = V x r(r, r') x V'. (10b) 

Thus, from the knowledge of the classical Green's dyadics, we can calculate 
the vacuum energy or stress. 

We now introduce the appropriate partial wave decomposition for a cylinder, 
in terms of cylindrical coordinates (r, d,z), a slight modification of that given 
for a conducting cylindrical shell [14] ^ : 



r'(r,r»= E / 1^ (Vxz)/^(r;A;,a;)x..fc(^,^) 

+ —V X (V X z)gmir;k,uj)xmk{d,z) \, (11a) 

LOE J 

dk ( 

*(r,r';u;)= J2 / t^UV x z)^^(r; /c, u;)Xmfe(^, ^) 

- —V X (V X z)fm{r;k,u;)xmk{0,z) I, (lib) 

LO/J, J 

where the cylindrical harmonics are 

X(^,^) = ^e^-^e^^ (12) 



and the dependence of fm etc. on r' is implicit. Notice that these are vectors in 
the second tensor index. Because of the presence of these harmonics we have 



^ It might be thought that we could immediately use the general waveguide decom- 
position of modes into those of TE and TM type, for example as given in Ref. [27]. 
However, this is here impossible because the TE and TM modes do not separate. 
See Ref. [28]. 
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r dr ' 

V X ( V X z) ^ Y\k- zdrt 

or r 



AT, 



in terms of the cylinder operator 

dm. — 



1 d d 



r dr dr r^ 

Now if we use the Maxwell equation (5b) we conclude ^ 



(13b) 



(14) 



^2\^ _ ,.2. 



(dm - k )fm^-^ l^fm- 



(15a) 
(15b) 



Prom the other Maxwell equation (5a), we deduce (we now make the second, 
previously suppressed, position arguments explicit; the prime on the differen- 
tial operator signifies action on the second primed argument) 



dmVjmir; r', 9\ z') = "^M'^Uir - r')xU{e', z'), 

e r 

dm-Dmgmir; r', 9', z') = -\ujJ^'*U{r - r')Xmk{0' , z'), 

where the Bessel operator appears, 

T^m — dm + ^ 1 ^ — i^'^SfJ' — k^- 
Now we separate variables in the second argument, ^ 



(16a) 
(16b) 

(17) 



L(r, r') = M'*Fm{r, r'; k, cu) + -M'*Fm{r, r'; k, uj) 



UJ 



(r, r') 



M'*Gm{r, r'; k, uj) - iM'*G„(r, r'; k, uj) 



X*mk{0\z'), (18a) 



where we have introduced the two scalar Green's functions F^, G^, which 
satisfy 



^ The ambiguity in solving for these equations is absorbed in the definition of 

subsequent constants of integration. 

^ Note that here and in the fohowing there are sUght changes in notation, and 
numerous corrections, to the treatment sketched in Ref. [1]. 
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dmVmFmir, t') = ^-5(r - r'), (19a) 
£ r 

dmVmGmir, r')^u;''^5(r - r'), (19b) 

while and are annihilated by the operator d^V^, 

dmVmF{r, r') = dmV^G{r, r') = 0. (20) 

The Green's dyadics have now the form: 

°° r°° dk ( / d - P\ 

+ -MAT" f-^^^) + A/-Ar"4-G„(r,r') 

+ —MM'*Grn{Ty)\xrnk{e,z)x*rnk{0\z'), (21a) 



'r»= E / 7^ --MAA'*G^(r,r')-iMM'*G^(r,r') 
- —AfM'*Fm{r,r') - -^ArAf'*Frn{r,r')\ 

CU/I iV^/I J 

xxmk{e,z)xU{e\z'). (21b) 



In the following, we will apply these equations to a dielectric-diamagnetic 
cylinder of radius a, where the interior of the cylinder is characterized by 

a permittivity e and permeabihty fi, while the outside is vacuum, so e = 
fi = 1 there. Let us consider the case that the source point is outside, r' > 
a. If the field point is also outside, r, r' > a, the scalar Green's functions 
F^, G'^, F', G' that make up the above Green's dyadics (we designate with 
primes the outside scalar Green's functions or constants) obey the differential 
equations (19a), (19b), and (20) with e = = 1. To solve these fourth-order 
differential equations we introduce auxiliary Green's functions ^'^ ■* (r, r') 
and ^^ '•'^'•'(r, r'), satisfying (m 7^ 0) 



dmg^^'''\r,r') = ^S{r-r'), (22a) 
dmg^^'^'\r,r') = Q, (22b) 

which therefore have the general form 
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-, / \ |m| 

er'-'(r,0^«-<«>(/)^-^(r^) . (23a) 

6r(a''(r,r')=af<a)(r')^, (23b) 

where r<(r>) is the lesser (greater) of r, r' and we discarded a possible r'™' 
term because we seek a solution which vanishes at infinity. Thus F^, G'^, F' 
and G' satisfy the second-order differential equations 



V^F'^^uj'gC. V^G'^^u^g^i (24a) 
P^F; = VrnG'^ = a;2g,G' (24b) 

Now, from (24a) and the first identity in (17) we learn that (A'^ = cu^ — k"^) 

K - = K{r')HUX'r) - y,jJM'r<)HUX'r^), (25) 

while G'^ obeys a similar expression with the replacement F ^ G. Similarly, 
fi-om (24b) 

FL-Y^i-<{r')Hm{X'r), (26) 

and for G'^ replace F ^ G. Here, to have the appropriate outgoing-wave 
boundary condition at infinity, we have used Hm{\'r) — H^\X'r). 

The dependence of the constants on the second variable r' can be deduced 
by noticing that, naturally, the Green's dyadics have to satisfy Maxwell's 
equations in their second variable. Thus, by imposing the Helmholtz equations 
in the second variable together with the boundary conditions at r' = oo, it is 
easy to see that 



<ir')^a'^-^^+b'ZHUX'r'), (27a) 
= + B'^H^X't'), (27b) 

and with similar relations for a'^{r'), A'^{r'),a'^{r'), and so on. Then, the 
outside Green's functions have the form 



KM- 



u 
+ 



/\m\ 



-\m\ 



id 



1 



— 

A'^ 2\m\ \ r> 



|m| 



A/F 



2 

(28) 
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while has the same form with the constants a'^, b'^, A'J[, B'^ replaced by 
a'^, b'^, A'^, B'^, respectively. The homogeneous differential equations have 
solutions 



rl\m\ 



(29) 



while in G'^ we replace a a"^, etc. 



When the source point is outside and the field point is inside, all the Green's 
functions satisfy the homogeneous equations (20) with e, /i ^ 1, and then, 
following the above scheme we have that 



and 



since now r can be 0. Also VmFm = ^"^Qm ^-^d 



(30) 
(31) 



(32) 



Gm,Fm and Gm have the same form, and the constants a^(r'), A^(r'), etc. 
follow the pattern in (27a) and (27b). Now, we may write for r < a,r' > a 



UJ 



— + b^H (X'r') 



(33) 

and similarly for Gm, F^, Gm, with the corresponding change of constants. In 
all of the above, the outside and inside forms of A are given by 



(34) 



The various constants are to be determined, as far as possible, by the boundary 
conditions at r = a. The boundary conditions at the surface of the dielectric 
cylinder are the continuity of tangential components of the electric field, of 
the normal component of the electric displacement, of the normal component 
of the magnetic induction, and of the tangential components of the magnetic 
field (we assume that there are no surface charges or currents) : 

is continuous, eEn is continuous, 
Hj is continuous, /iHn is continuous. (35) 

These conditions are redundant, but we will impose all of them as a check of 
consistency. In terms of the Green's dyadics, the conditions read 
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e ' v 

r-sT' 

e • * 

z • $ 



r=a+ 



r=a— 
r=a+ 



r=a— 
r=a+ 



r=a— 

r=a+ 



r=a— 
r=a+ 



r=a— 
r=a+ 



0, 
0, 
0, 
0, 
0, 
0. 



(36a) 
(36b) 
(36c) 
(36d) 
(36e) 
(36f) 



We can also impose the Helmholtz equations (9a) and (9b). From those we 
learn that the coefficients of terms with powers of r are related in the following 
way 



a'P + a'G ^ 0, (37a) 
6'^-(sgnm)-6'^ = 0, (37b) 

6'^- (sgnm)-6'^ = 0, (37c) 

for the Green's dyadics outside the cyhnder and equivalent expressions for the 
inside (no primes) 



ia^-a^ = 0, (38a) 
A* 

6^ + (sgnm)--6^ = 0, (38b) 

\P + (sgnm) - -6^ = 0, (38c) 

where we have introduced the abbreviations for any constant K 

^K"" - (sgnm) -K^", ^ - (sgnm) "^K^ , (39) 

oj k 

and the same for K'^ and K'^ (the outside). Then, from the boundary condi- 
tions we can solve for the remaining constants. Notice that, due to the tensorial 
character of the Green's dyadics, each of the above six boundary conditions 
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(36a). (36b), (36c), (36d), (36c), (36f) arc in fact three equations correspond- 
ing to the three prime coordinates. From (36a) we get the following three 
equations: 



- e\aJ'^{\a)B^ Jm{\a)Bg = -\'aH'^{\'a)B'^ 

CUE 

H,,{X'a)B'^ + ^-JM'a), 40a 

-s\m\XaJ^{Xa)A^ H XaJ^{Xa)A^ H —Jm{Xa)A]^ 

UJ UJ £ 

-'^UXa)Ai = -\m\X'aHi{X'a)A': + ^X'aH'JX'a)AZ 

+ ^HM'a)A'^ - ^^Hrr.{X'a)A'i, (40b) 
eXaf^{Xa)B^ + —UXa)Bg = X'aH^iX'a)B'^ 

UJE 

- ^|A'aj;(A'a) + ^H^{X'a)B'^. (40c) 



The three equations following from (36b) are: 



\m\A, 



G 



B 



G 




2 r 



,QHm{X'a) lJ^ TT JmjX'a) 

A'2 2i Jm{Xa) 



B 

Jm{Xa) 
Hm{X'a 



Jm{Xa) 



UJ 



iG 



m\A'^ - m-A^ 



(41a) 
(41b) 
(41c) 



and those coming from (36c) are: 
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me^Jm{\a)B^ + -XaJ'^{\a)Bg = mHm{\'a)B'^ + -X'aH'^{X'a)B'^ 

00 LO 



+ 



mk 



m\m\k o ^ , , . , p k"^ 
m\m\k 



(42a) 



\m\ 



\aJ'^{\a)A 



LO 



\aJ'J\a)Ai^m^E^{}!a)J^l - ^=^//„(A'a)A 

LO 



IF 



^X'aH'JX'a)A'^ + '^X'aH'JX'a)At 

LO LO 



(42b) 



7 2 
e^mJm{Xa)B^ + -XaJ'JXa)Bg = mH^{X'a)B'i - J^(A'a) 

+ h'aH:^{X'a)B'l (42c) 

LU 



Prom the set of equations involving the magnetic part, we find that (36d) 
gives us 



ek 



ixm\m\k 



IJ,mJm{Xa)B^-\ XaJ'„X^a)B^ = mHm{X'a)B'^ 

LO 

-m^|j^(A'a) + ^A'a//;(A'a)sf, 

lim'JU\a)AZ - —XaJ^{Xa)A^ 



(43a) 



Jm{Xa)A, 



£|m|A;^ 



LO^ 



XaJ^{Xa)A: 



m|m|A; 



LO 

i'G _2; 



Hm{X'a)A'^ - m^Hm{X'a)A 
m\k'^ f 



LO 

^X'aHU{X'a)A': + 

LO 

ek 



iG 

m 



LO^ 



X'aH'^{X'a)A'aA3h) 



—IJ,mJm{Xa)B^ XaJ'^{Xa)B^ = —mHjn{X'a)B. 



jiG 



LO 



+ ^|A'aj;(A'a) - ^AW;(A'a)s;^. 



(43c) 



By imposing (36e) we get the conditions 
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XaJ'^{Xa)Bg + MXa)B^ = X! aH'^{X! a)B'^ 

+ —H^{X'a)B'i - ^^A'aj;(A'a), (44a) 

I TJli I /c ~ TTli^ 

-XaJ'^{Xa)Ag + mXaJ'^{Xa)Ag + Jm{>^a)A^ 



III, \ / III, III, \ / III, 

'^\'^\'^'-UXa)Ai = -^A'ai7;(A'a)A:^ + mX' aH'^{X: a) 



+ ^HM'a)AZ - ^^^^HM'a)A'l (44b) 

Ul UJ 

XaJ'^{Xa)Bi+'^Jr^{Xa)Bl = >^aH'^{X'a)B'i 

UJ/J, 

-^l:Uya) + ^HUX'a)B'^. (44c) 



And finally (36f) gives us 



p k \m\ p 




Hm{X'a) 
Jm{Xa) ™ 

H^{X'a) 



Jm{Xa) 



2 r 



k |m| 
u! m 

2 



A, 



IF 



B, 



,p Hjn{X'a) _ UJ^ TT JmjX'a) 

Jm{Xa) X'^2i Jm{Xa) 



(45a) 
(45b) 
(45c) 



By combining these equations we find the remaining constants, but the equa- 
tions are not all independent. First, from (41b), (45b), (40b) and (42b) we 
learn that the coefficients of terms involving Bessel functions and r'"'"^' cancel 
among themselves in a way such that the ones from the outside do not mix 
with those from the inside: 



4^ = 4^? = n 



(46a) 
(46b) 



The same can be found if we use (44b) and (43b) instead of (40b) and (42b) . 

Next we determine the coefficients of functions involving just Bessel functions. 
Prom (44c) and (42c) we find using (45c) and (41c) that 
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BZ^--{l-et^)^U\a)H^{^a)B^^, (47a) 
all in terms of 

^--—eXy'E ^^^^ 
found by subtracting ^ times equation (42c) from (44c) and using (43c) . The 
denominators occurring here are^ 



5= (1 - ef^r^^:^J'Ma)HUX'a) - DD, (49a) 

D = sX'aJ^(Xa)Hm(X'a) - XaH'^(X' a) Jm{Xa) , (49b) 
D^liX'aJ'^(Xa)Hm(X'a) - XaH'^(X'a)Jm(Xa). (49c) 

The second set of constants is found using (43a), (40a), (45a) and (41a): 



B^^-^^il-en)^UXa)HUX'a)BZ (50a) 

= - (^)' J(l-^/^)^^^(Aa)5«, (50b) 

G_ o;^ 71 JmjX'a) 1 JrniXa) Q 

"""^ X'^2iH^{X'a)^\X') eH^yaf"^' ^^^"^^ 

in terms of 

coming from (50b) and (42a) . 

It might be thought that m = is a special case, and indeed 



m\ 

^In^, (52) 



1 , T< 



2\m\ V r> / 2 r 



^ (43c) is the same equation as (40c), which can easily be seen by using (45c). 
^ The denominator structure appearing in E is precisely that given by Stratton 
[28], and is the basis for the calculation given, for example in Ref. [10]. It will be 
employed in an independent rederivation of the Casimir energy for a dilute dielectric 
cylinder [29]. 

^ By using (45a) it can be seen that this equation is the same as (44a). 
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but just as the latter is correctly interpreted as the limit as |m| ^ 0, so the 
coefficients in the Green's functions turn out to be just the m = limits for 
those given above, so the m = case is properly incorporated. 

It is now easy to check that, as a result of the conditions (37a), (37b), (37c), 
(38a), (38b), (38c), (46a), and (46b), the terms in the Green's functions that 
involve powers of r or r' do not contribute to the electric or magnetic fields. 
As we might have anticipated, only the pure Bessel function terms contribute. 
(This observation was not noted in Ref. [14].) 



3 Stress on the Cylinder 



We are now in a position to calculate the pressure on the surface of the cylinder 
from the radial-radial component of the stress tensor 



P=(T,,)(a-)-(r,,)(a+) 



(53) 



where 



T — - 



e{El + El - El) + ii{Hl + Hi - Hi) 



(54) 



As a result of the boundary conditions (35), the pressure on the cylindrical 
walls are given by the expectation value of the squares of field components 
just outside the cylinder, therefore 



-'-^\El + El + ^ 



(55) 



r=a+ 



These expectation values are given by (10a), (10b), where the latter may also 
be written as 

(H(r) H(r')) = -— *(r, r') x V'. (56) 

UJjJl, 



It is quite straightforward to write the vacuum expectation values of the fields 
occurring here outside the cylinder in terms of the Green's functions. 
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m=—oo 



mk d d^-k'^ ~ k'^ d d 

7^ 2 — ^rnir, r + —^^G^{r, r) 

ujr or uj'^ Lo^ or or 



m=—oo 



mk d dm-k^ p, , rr^ , , 



1 ^TTl ^_ J-oo ZTT CU^ 



According to (56) the magnetic field expectation values can be written as 
follows, 



{Hr{r)Hr{r'))^-^^ f H ^\ - ^{d'^ - k^)G'Jry) 
27r uj ^t^^ J-^ 27r [ ujrr' 

"9 ~, , k'^ d d ^, . , 

(58a) 

fj \ oo foo c\k ( 1 r? r? 

mk d ~, , ,, m^/c^ ^, , ,, 

{H,ir)H,ir'))^--- ^ / ^-^m«(r,r'). (58c) 

ZTT a; J-oo ZTT a; 



m=— oo 



When these vacuum expectation values are substituted into the stress expres- 
sion (55), and the property of d^ exploited, 

ci^r±H^O, dmJm{Xr)^-X^Jm{Xr), (59) 

(of course, the later formula holds for Hm as well and the same for d'^ acting 
on the primed coordinate), we obtain the pressure on the cylinder as 
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47ri ^±^^7-00 27r i- 



°° do; /-"^ dkX^ 



+ i/;(a;')4(^)^m(^') 



(2e + 2)(l-£/i) 



+ 



A2 



(1 - enf + xW 



m 



X'- 



fc^ + — + A' 



/2 



- E'Jx')Jl{x)E^{x')ix}!^x\u?e + 



m 



X' 



(60) 



where x — Xa, x' — X'a and the last bracket indicates that the expression 
there is similar to the one for the electric part by switching e and showing 
manifest symmetry between the electric and magnetic parts. 

In order to simplify this expression, we make an Euclidean rotation [30], 

a; — > iC A — > ire, (61) 
so that the Bessel functions are replaced by the modified Bessel functions. 



Jm{x)Hm{x') —Im{y)Km{y'), 
TTl 



(62) 



where y — Ka and y' — n'a. Then (60) becomes 



£-1 



E 



167r3a4^_^^-oo 



- K'M)ll{y)KUy') 



dCa dkaii^K'^{y')I^{y)I'^{y)y{k'^a'^ - Ca^i^) 
w?'k'^a?C,'^a? ^ 



y'^e 



2(£ + l)(l-£/i) 



+ 



y^ ( / 2 2 C^a^\ , ^,,2 

K Qj 



y' \y'^ 



Uy)Uy)Kliy')y 



m 



::^(^V/.-CV) + y'V 



+ y'- 



+ {e^ fj,), 



(63) 



where 
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5 = lliy)Kl{y'){l - e^f + AA, (64a) 

y y 

A = ey'Ii{y)KUy') - yK'^{y')Im{y) (64b) 

K = ,,y'I'^{y)KM) - yK'Jy')I^{y) (64c) 



This result reduces to the well-known expression for the Casimir pressure when 
the speed of light is the same inside and outside the cylinder, that is, when 
eyu — \. Then, it is easy to see that the denominator reduces to 



S = AA 



4£ 



/12 



1 - ey\i^Kj\ 



(65) 



where ^ = {e — !)/(£ + 1). In the numerator introduce polar coordinates. 



= /c^a^ + C^a^, ka — y sin 9, (a — ycos9, 



(66) 



and carry out the trivial integral over 6. The result is 



1 POO ^ A 

■^i <iyy\Z^^H^-eMKMf), (67) 



which is exactly the finite result derived in Ref. [10], and analyzed in a number 
of papers [11,12,13]. 



4 Bulk Casimir Stress 



The expression derived above, (63), is incomplete. It contains an unobservable 
"bulk" energy contribution, which the formalism would give if either medium, 
that of the interior with dielectric constant e and permeability /x, or that of 
the exterior with dielectric constant and permeability unity, fills all the space 
[7]. The corresponding stresses are computed from the free Green's functions 
which satisfy (19a) and (19b), therefore 



2 m V r> / 2i 



(68) 

where < r,r' < oo. Notice that in this case, both F^^ and are zero and 
the Green's dyadics are given by 
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+ -^MAf'*Gt\Ty)\xnrk{0,z)x*^^{e\z'), (69a) 

°° roo Hi- ( i 



ArM'*Fj^\r, r') \xmk{0, z)x'mk{0\ z'). (69b) 



It should be noticed that such Green's dyadics do not satisfy the appropriate 
boundary conditions, and therefore we cannot use (55), but rather one must 

compute the interior and exterior stresses individually by using (54). Because 
the two scalar Green's functions differ only by a factor of n/s in this case, for 
the electric part the inside stress tensor is 



da; 



dA; 1 



+ \-d' 



m 



rr 



-oo 27r Oi)^£ 



d d 



(70) 



r=r—a— 



while the outside bulk stress is given by the same expression with A — > A' = 
uj^ — and e = ^ = 1. When we substitute the appropriate interior and 
exterior Green's functions given in (68), and perform the Euclidean rotation, 
Lu —>■ i(, we find a rather simple formula for the bulk contribution to the 
pressure 



P^ = r/o)(a-)-T,(°)(a+) 

h 



167r3a4 



E 



dCa dka\ y'liiy)K'M - (r + m')Im{y)KM 



- y"Uy')Kiy') + (r + m')Uy')KM) 



(71) 



This term must be subtracted from the pressure given in (63). Note that this 
term is the direct analog of that found in the case of a dielectric sphere in 
Ref. [6]. Note also that P*" = in the special case eii — l. 

In the following, we are going to be interested in dilute dielectric media, where 
H = I and |£ — 1| -C 1. We easily find that when the integrand in (71) is 
expanded in powers of (e — 1) the leading terms yield 
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87r2a4 



+ 




(72) 



where in the last form we have introduced polar coordinates as in (66) and 
performed the angular integral. 



5 Dilute Dielectric Cylinder 



We now turn to the case of a dilute dielectric medium filling the cylinder, that 
is, set /i = 1 and consider e —1 as small. Wc can then expand the integrand in 
(63) in powers of {e — 1) and, because the expression is already proportional 
to that factor, we need only expand the integrand to first order. Let us write 
it as 



where we have noted that the (e — 1)^ in S (64a) can be dropped. Expanding 
the numerator and denominator according to 



N = Ar(o) + {6- l)iV(^) + . . . , AA = 1 + (£ - l)A(i) + . . . , (74) 




(73) 



we can write 



P 



{e-l)h 




where 
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7V(0) = -(^V-CV)X;(y)7;(y) 



m 



Km{y)Im{y)-, 



2^2 



1 + 



m 



(fcV-CV)i^-(y)4(i/) 



2 

_cv 

2 

+ CV 



1 + 



y 



4 



y ( 1 + ^ ) + 5-(fc^a^ - Ca') - -,m'k'a' 



+ 
+ 



m \ m 

^2 J y 

/•2 2 I S " / 7 2 2 a2 2\ 

y 

i2\ 



(76a) 



Km{y)Kiy)Ciy), 



(76b) 



+ CV|1 + 



m 



r 



lm{y)Km{y)- 



(76c) 



When we introduce polar coordinates as in (66) and perform the trivial angular 
integrals, the straightforward reduction of (75) is 



h 



[e-l) Yl / dyWKm{y)Ir, 



iy) 



^K'^iy)Uy)Uy) 



+ K'Jy)Ciy)KUy) + K'^{y)C{y)l 



m 



Kiymyy-: 1 + + Kiymyy,; 1 + 



m 



+ Kl{y)i'Mim{y) (i + ^) 



m 



m 



KL{y)C{yr- ( 1 - ^ 



(77) 
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The leading term in the pressure, 



^ / ^yy'Km{y)Uy). (78) 



can also be obtained from (63) by setting e — ^ — 1 everywhere in the 
integrand, and the denominator S is then unity. This is also exactly what is 
obtained to leading order 0[(e — 1)^] from the bulk stress (72). Thus the total 
stress vanishes in leading order: 

P«-P^ = 0[(£-1)'], (79) 



which is consistent with the interpretation of the Casimir energy as arising 
from the pairwise interaction of dilutely distributed molecules. 



6 Evaluation of the (e — 1)^ term 



We now turn to the considerably more complex evaluation of the {e — lY term 
in (77). 



6.1 Summation method 



As a first approach to evaluating this second-order term, we first carry out 
the sum on m by use of the addition theorem 



MkP)^ Y: e'-('^-*')X^(M/m(V), P>P', (80) 



where P = y p"^ + p'"^ — 2pp' cos(0 — 0'). Then by squaring this addition theo- 
rem and applying suitable differential operators, in the singular limit p' — > p 
we obtain the following formal results: 
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5: Ki{kp)ii{kp) 



m=—oo 
oo 



2n 



27r 



\2 2 
^ COS 



m=:^cx) -^0 27r 2 

oo 



2^ d0 
27r 



^ Kl(kp)Ukp)Ukp) 

m=—oo 
oo 

Y: m'Kl{kp)Ukp)Ukp) 

=— oo 

oo 



2^ d02; 2 

^2 2' 



m=— oo 
oo 



2lT 



Ztt I 



E ^'Kl{kp)ll{kp)=r ^\kI{z)^o.% 
„^=^oo Jo 27r 4 2 

£ Ukp)lLikp)K'^{kp)^ r ^K{z)sm^ 
.^r^ -fo 27r 2 



(81a) 
(81b) 

2 

2J ' 

(81c) 
(81d) 

(81e) 

(81f) 

(81g) 



X 



Ko{z) sin^ - 



K'oiz) 



00 2tt 

E I'l{kp)K'^{kp)^ 



d0 
2^ 



Kq[z) sm 

z z 



(81h) 



(81i) 



Here z = 2A;psin |, and we recognize that in this singular hmit (which omits 
delta functions, i.e., contact terms) terms with and interchanged in 
the sum have the same values. (For further discussion of this, see Ref. [29].) 

When we put this all together, we obtain the following expression for the 
pressure at second order: 



P(2) = 



40967r2c 



r dzz' r ^<t>\ K{z)-rKl{z){l-A/z^) 
r^a^Jo Jo 27r\ sin^ (J)/2 



^ ^ {1 - 8/ z')K',iz) -2(1 + 3/ z')K',\z) Kl{z) 



is -I) 



153607r2a4 Jo 2-k 



2tt 



sinV/2 
5 



z* sin^ (j)/2^ 
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sin 0/2 sinV/2 sin^0/2 



(82) 



Of course, the (f) integrals in (82) are divergent. However, we will regulate them 
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by continuing from the region where the integrals converge: 




(83) 



which is vahd for Res > —1. We will take the right side of (83) to define 
the angular integral for negative s. Then wc sec that those integrals vanish 
when s = —2n where n is a positive integer. Thus, this analytic continuation 
procedure says that the result (82) is zero. As for the bulk term, the addition 
theorem (80) implies that the y integral in the second term in (72) reduces to 



This argument is exactly that given in Ref. [2] to show that the Casimir energy 
of a dilute dielectric-diamagnetic cylinder with sfi = 1 vanishes. However, it 
is not very convincing, because it seems to show no relevance of cancellations 
between various terms in the expressions for the pressure. That relevance will 
be established in the method which follows. 



6.2 Numerical analysis 

We now turn to a detailed numerical treatment of the second-order terms in 
(77) and (72). It is based on use of the uniform asymptotic or Debye expansions 
for the Bessel functions, m » 1: 




(84) 




(85d) 



(85b) 



(85a) 



(85c) 



where y = mz and t = + z^. (The value of r] is irrelevant here.) The 

polynomials in t appearing here are generated by 
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uo{t) = 1, vo{t) = 1, (86a) 

Uk{t) = - t'')ui_,{t) + 1 f dx{l- 5x'')uk-i{x), (86b) 
2 8 JO 

Vkit) = nfe(i) + t{t^ - 1) (^^ife-i(i) + K-i W) ■ (86c) 
Now suppose we write the second-order expression for the pressure as 

^=W^i;;i d„vw. (87) 

where the exphcit form for gmiv) can be immediately read off from (77) and 
(72), and the prime on the summation sign means that the m = term is 
counted with half weight. We have recognized that the summand is even in 
m. Let us subtract and add the first five terms in the uniform asymptotic 
expansion for g^^, m ^ 1: 

sM-^i^M^), (88) 

where z — y/m and 



fi{z 

f2{z 

h{z 

Mz 



4^(1 + ^2)3' 

-8 + 8^2 + ^4 

16 - 84;^^ ^ 34^4 _ iQ^Q 



5z' 



16^(1 + ^2)6 

-64 + 1024^2 _ 1854^4 ^ ^Q^^t 

64z(l + ^2)13/2 

64 - 2416^2 + 11808^4 _ 15595^6 



(89a) 
(89b) 
(89c) 
(89d) 



6856^* - 555z^° - 15z 



vl2 



64^(1 + ^2)5 



(89e) 



We note that when these functions are inserted into (87) in place of Qm, the 
first three fk give divergent integrals, logarithmically so for /i and fs, and 
linearly divergent for /2. We also note the crucial fact that 



r dzz^u{z)^o, 

Jo 



(90) 



which means that C{1), which would indicate an unremovable divergence, does 
not occur in the summation over m. This is the content of the proof that the 
Casimir energy for a dilute dielectric cylinder is finite in this order, given by 
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Bordag and Pirozhcnko [24]. We also note that when the divergent part is 
removed from the /2 integration we again get zero, 

f d.(.V.W-i)=o. (91) 

The suggestion is that this term may be simply omitted as a contact term. 
(But see below.) 

However, the two logarithmically divergent terms, corresponding to /i and /s, 
give finite contributions, because they are multiplied by formally zero values of 
the Riemann zeta function. The first one may be regulated by a small change 
in the power: 

hm Y: m'-' / dzz'-^Mz) = ]im C{-2 + s)- = (92a) 
The term gives similarly 

lim E dzz'-^Mz) = C'(0) {--) = -ln27r. (92b) 



Although it would appear that a finite term would emerge from f^, that term 
vanishes because remarkably 

/■oo 

/ dzzUiizfJz)=0. (93) 
The /5 term is completely finite: 

] POO 1 07^2 

dzz^Mz) = ^. (94) 
^^m^ Jo ^ ' 7680 ^ ' 



Following the above prescription, we arrive at the following entirely finite 
expression for the pressure on the cylinder: 



327r2a4 \ 

00 

m=l 



C(3) 5 , ^ 197r2 
+ 7^ In 27r + 



167r2 32 



7680 



9m{y) 



2m^ ^ m 



+ 



/■oo 

/ dy\ 
Jo 



9o{y) - - \h{.y) 



(95) 



Here in qq we have subtracted a linearly divergent term, which when combined 
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with that removed in (91) gives 



100 

^ E 7 d?/- (96) 



We regard this, rather cavaherly, as a contact term, which we simply omit. In 
the next section we will give the correct treatment of this /2 term. All that 
remains is to do the integrals numerically. We do so for m from through 4, 
after which we use the next nonzero term in the uniform asymptotic expansion, 



because, again, the integral over /e vanishes. 

When all the above is included, to 6 decimal places, we obtain 



is - I? 

P = ^ ^(-0.007612 + 0.287168 + 0.024417 - 0.002371 - 0.000012 

327r^a* 

- 0.301590) = 0.000000, (98) 

where the successive terms come from (92a), (92b), (94), the numerical integral 
over the first 4 subtracted g^s (m > 0), the remainder (97), and the numerical 
integral over the subtracted go, respectively. This constitutes a convincing 
demonstration of the vanishing of the Casimir pressure in this case. It is similar 
to the numerical demonstration [10] of the seemingly coincidental vanishing 
of the Casimir energy for a dilute dielectric-diamagnetic cylinder, obtained by 
expanding (67) to order 



6.3 Exponential regulator 



Although the calculation in the previous subsection is quite standard, and 
undoubtedly correct, the reader might rightly object that zeta-function regu- 
lation has been employed, and infinite terms simply omitted. Therefore, and 
to make contact with known results, let us insert a regulator to make all the 
sums and integrals completely finite. It would be best, as in Ref. [14], to insert 
such a regulator before rotating the frequency in the complex plane. However, 
this is much more complicated here than in that reference; and because the ex- 
pressions here are formally much more divergent, the regulator adopted there 
appears insufficient. It will suffice for the present purposes to simply insert by 
hand an exponential regulator into the expression (87): 

_ 1 ';2 00 „oo 

P...= ^-^j:;i d„Vto)e-», (99) 
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where S — > 0+ at the end of the calculation. Then it is easy to repeat the 
calculation of the previous subsection. One has only to carry out the sum 

oo -1 

E e-^"^^ = (100) 

m=l ^ ^ 

Then the /i term, instead of (92a), is 

r ^ 4. . ^ 1 137r C(3) 

i ^^^^^^^)d^i^=32^-li^- ^'''^ 

The /2 term has no finite part: 

/■oo HI 1 

d^zV2(^)^^^ — = (102) 

io ^ ' d6ze^' - 1 16(5' ^ ^ 

where the reader should note that no ad hoc subtraction as in (91) has been 
employed. The evaluation of (102) uses the fact that 

/■oo 

/ d^^V2(-2) = 0. (103) 
Jo 

The /s term is, instead of (92b), 

/; d-V3(.) + ie-) ^ + |b.2.. (104) 

Here we have subtracted a term from the m = contribution: 

foo r 1 1 1 

dyy' \go{y) - -fs{y)\ e"^^ = -0.301590 + — . (105) 

The divergent term here cancels that in (102), and the finite part is the value 
of the last integral in (95). Thus wc recover exactly the same numerical result 
(98) found in the previous subsection, plus two divergent terms 

_ (^-1)^ A^TT^ 3157r\ 

~ 1, 32^ ~ 8192^ ) ■ ^ ' 

The form of the divergences is exactly as expected [24,23]. In particular, there 
is no 1/5^ divergence, because of the identity (103). 



6.4 Interpretation of divergences 



In the previous section we computed divergent contributions to the Casimir 
pressure for a dilute cylinder. For simplicity, we chose an exponential regulator 
with a small dimensionless parameter 5 — > 0+. How do we interpret these 
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terms? It is perhaps easiest to imagine that S as given in terms of a proper- 
time cutoff, 5 = T / a, T ^ 0+. Then if we consider the energy, rather than the 
pressure, the divergent terms have the form 

^ aL LI , ^ 

Ediv = es— + ei . 107 

Here L is the (large) length of the cylinder. Thus, the leading divergence 
corresponds to an energy term proportional to the surface of the cylinder, and 
it therefore appears sensible to absorb it into a renormalized surface energy 
which enters into a phcnomenological description of the material system. (Such 
arguments are familiar, dating back to [6], and recently vigorously revived in 
[31].) The 1/r divergence is more problematic. It is proportional to the ratio 
of the length to the diameter of the cylinder, so it seems hkely that this 
would be interpretablc as an energy term referring to the shape of the body. 
If one could compute the Casimir energy of an extremely elongated ellipsoid, 
we would expect an energy term proportional to the ratio of curvatures. (Of 
course, a cylinder has zero curvature.) This appears to be exactly of the form 
of a surface integral [16] 

j dSKiK2, (108) 

in terms of the principal curvatures Ki, i = 1,2. Siich terms are well known 
not to contribute to the observable energy. Had a divergent term proportional 
to appeared in the pressure, it would have implied a divergent energy of 
the form 62 (In a) /r^, which would have been impossible to remove. (For the 
dielectric sphere the situation is simpler, in that divergences are all associated 
with positive powers of the sphere's radius [8].) 

In any case, although the structure of the divergences is universal, the coeffi- 
cients of those divergences depend in detail upon the particular regularization 
scheme adopted. In contrast, the term proportional to {e — l)^/a^ is unique. 
Thus, of course, it could not have been any other than that zero value given 
by the van der Waals calculations [10,22,21]. 

The nature of divergences in such Casimir calculations is still under active 
study [19,20,1]. The universality of the finite Casimir term makes it hard not 
to think it has some real significance. As an example of how subtle interpreta- 
tion of divergences can be, we recall that it has now been proved that the total 
Casimir energy for electromagnetic modes interior and exterior to an arbitrar- 
ily shaped smooth infinitesimally thin closed perfectly conducting surface is 
finite [32] . This is hard to reconcile with the existence of local divergences in 
the energy density near the surface proportional to (ki — ^2)^. Presumably, 
these divergences belong to the surface itself and have nothing to do with 
the global Casimir energy [20,1]. But the open questions are profound and 
challenging. 
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7 Conclusions 

Since the beginning of the subject, the identity of the Casimir force with van 
der Waals forces between individual molecules has been evident [33,34]. It is 
essentially just a change of perspective from action at a distance to local field 
fluctuations. So it was no surprise that the retarded dispersion force between 
molecules, the Casimir-Polder force, could be derived from the Lifshitz force 
between parallel dielectric surfaces [35,26]. However, the identity is not really 
that trivial, because both the van der Waals and the Casimir energies contain 
divergent contributions. This is particularly crucial when one is considering 
the self-energy of a single body rather than the energy of interaction of distinct 
bodies. Thus it was nontrivial when it was proved that the Casimir energy of 
a dilute dielectric sphere [8] coincided with that obtained by summing the van 
der Waals energies of the constituent molecules [9]. 

When it was shortly thereafter discovered that the sum of van der Waals forces 
vanished for a dielectric cylinder [21,10] it was universally believed that the 
corresponding Casimir energy, in the dilute approximation, must also vanish. 
This result was verified by a perturbative calculation [23]. Proving this by a 
full Casimir calculation turned out to be extraordinarily difficult. This paper 
is the result of a five- year-long effort. It should dispell any lingering doubts 
about the meaning of the Casimir force. The importance of this finding is 
impossible to evaluate at this point; a zero value suggests some underlying 
symmetry, which is certainly far from apparent. It probably has technological 
implications, for example in the physics of nanotubes, which will be explored 
in a subsequent publication. 
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